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ABSTRACT 

Recent work has raised the exciting possibiUty that r-modes (Rossby waves) in rotating neutron star 
cores might be strong gravitational wave sources. We estimate the effect of a solid crust on their viscous 
damping rate and show that the dissipation rate in the viscous boundary layer between the oscillating 
fluid and the nearly static crust is > 10^ times higher than that from the shear throughout the interior. 
This increases the minimum frequency for the onset of the gravitational r-mode instability to at least 
500 Hz when the core temperature is less than 10"'^° K. It eliminates the conflict of the r-modc instability 
with the accretion-driven spin-up scenario for millisecond radio pulsars and makes it unlikely that the 
r-mode instability is active in accreting neutron stars. For newborn neutron stars, the formation of 
a solid crust shortly after birth affects their gravitational wave spin-down and hence detectability by 
ground-based interferometric gravitational wave detectors. 

Subject headings: dense matter - gravitation - stars: neutron - stars: rotation — stars: oscillations 



1. INTRODUCTION 



Th e recently discovered r-mod e instability ( Anderssorj 



1998; Friedman & Morsink 199S) might play a significant 



ing neutron stars ( 


Lindbloni ct al. 1998 


Owen et al. 1998; 


Andcrsson et al. : 


L999a; Bildsten 1998 


Andersson et al 


1999b). The unstable regime for the r-modes in the neu- 



tron star (NS) core depends on the competition between 
the g ravitational radiation e xcita tion and viscous dissipa- 
tion. Lindblom ct al. (1998D and Andersson et al. (1999a) 
computed the dissipation due to shear and bulk viscosities 
in normal fluids (not supcrfluids) and found that gravita- 
tional excitation exceeds vis cous damping in all but the 
most slowly rotating stars. Lindblom & Mendell (1999) 
showed that superfluid mutual friction is also not conipet- 
itive with gravitational radiation unless the superfluid en- 
trainment parameter assumes a very special value. Hence, 
on theoretical grounds, r-modes are expected to be unsta- 
ble over much of the parameter space occupied by newborn 
NSs, NSs in low-mass X-ray binaries (LMXBs), and mil- 
lisecond radio pulsars (MSPs) during spin-up. 

However, observations pose challenges to this theoret ' 



LMXBs (ran der Klis 1999), which would otherwise have 
been spun up to > 1000 Hz during their ^ 10^ yr lifetime 
( van Paradijs fc \yhitc 1995 ). However, for normal fluid 
cores. Levin (1999| ) showed that the temperature depen- 
dence of the shear viscosity makes this spin equilibrium 
thermally unstable, leading to a limit cycle behavior of 
rapid spin-down s after prolonged periods of spin-up. For 
superfluid cores. Brown fc Ushomirsky (1999 ) showed that 
such constant-amplitude spin-up is inconsistent with the 
quiescent luminosities of several known LMXBs. 

All but the hottest NSs have solid crusts that occupy 
the outer w km of the star. In this paper, we calculate 
a previously overlooked dissipation: the viscous boundary 
layer between the oscillating fluid core and static crust. 
The shear viscosity of NS matter is relatively small and 
does not affect the r-mode structure in the stellar interior. 
However, the transverse fluid motions of the r-modes are 
very large at the crust-core boundary and must, of course, 
go to zero relative velocity at the crust. The dissipation in 
the resulting viscous boundary layer (hereafter referred to 
as the VBL) substantially shortens (by t ypically 10^) the 
r-mo d c damping times. Previous work ([Lindblom et aT 



iPPi p.t.nr. Th. P^i.t.n^P nf twn 1 R rr.. raHin pnl.«r. 1998| [Liudblom ct al. 1999| ; |Andersson et al. 1999a) ) in 



means that rapidly rotating NSs are formed in spite of the 
r-mode instability (Andersson ct al. 1999b). While it is 



not clear that their current core temperatures place these 
MSPs inside the r-mode instability region, current the- 
ory says that they were certainl y unstable durin g spin -up. 
For accreting NSs in LMXBs, [Bildsten (1998D and [An] 
dersson et al. (1999t) conjectured that the gravitational 
radiation from an r-mode with a small constant ampli- 
tude could balance the accretion torque. This would pos- 
sibly explain the preponderance of « 300 Hz spins among 



eluded damping from the shear viscosity acting only on 
the gradient in the transverse velocity of the mode in the 
interior of the star (on the length scale of the stellar ra- 
dius, R). The shear in the VBL acts on a much shorter 
length scale, and is hence stronger. 

This new source of dissipation raises the minimum fre- 
quency for the r-mode instability in NSs with a crust to 
^ 500 Hz for T « 10^" K, and even higher frequencies for 
lower temperatures. It thus alleviates the conflict between 
the accretion-driven spin-up scenario for MSPs and the 
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r-mode instability a nd significantly alt ers the spin-down 
scenario proposed by Owen et al. (1998) for newborn NSs. 



2. THE VISCOUS BOUNDARY LAYER 

The transverse fluid motions of an r-mode cause a time- 
dependent "rubbing" of the core fluid against the other- 
wise CO- rotating crust spinning at f2 = 2'itVs. In previous 
works, the boundary condition applied at this location al- 
lowed the fluid to have large-amplitude transverse motion. 



To calculate the damping m ore accurately, we use the 
dissipation per unit area from Landau & Lifshitz (1959) 
and integrate over the surface area at the crust-core 
boundary 



—r = -\ -TT \^\ R sn 
dt J 2 \ 2 



(3) 



Negle cting viscosity is an eAcelleiil assumpLion fai 
from the crust core interface. — However, there can be no 
relative motion at the boundary for a viscous fluid, lead- 
ing to a VBL where the transverse velocity drops from a 
large value to zero (the "no-slip" condition). The VBL is 
mediated by the kinematic viscosity, in the fluid just 
beneath the base of the crust. The density there h as most 
recently be en estimated as pb ~ 1.5x10^"* g cm 



The fluid velocity is given hy v = (ar^i?)y,f^e'"* , and 
E ^ 1.64 X 10"^Af(arJi?)2/2 is the mode energy (pwcn 



et al. 199^ ). Integrating and using the mode frequency 
in the rotating frame, we find the damping time due to 
rubbing 



'^rub 



2E 
' dE/dt 



100 s 



1 kHz y/' 



et al. 199£). We use the values for v as found by 
fc Itoh (1979]) and fit by putler fc Lindblom (198 




X 10^ cm^ s 



Ti 



(1) 



In this fit, three different cases are handled with the pa- 
rameter /. When both neutrons and protons are nor- 
mal, the predominant scatterers are neutrons, and / — 
)^^^ , — If neutrons are superfluid and protons nor - 



(p/Pb 

mal, hhn visrnsity is mnrliatnrl hy nlectrnn-prnt.nn scat- 



(4) 

where p is the density at the crust-core boundary. 

So far, we have assumed that the large-scale relative mo- 
tion between the oscillating core and the crust is simply 
set by the toroidal motion of the r-mode. This is fine in the 
inviscid limit {v — 0). The transverse displacement at 
the crust-core boundary is then discontinuous, and the r- 
modes do not couple to the crust's torsional modes, which 

have ^ = ^±^im- ^'^^ 1^ — 0, r- modes can only couple to the 
crust via their 0{il^) radial motion. However, the VBL ex- 
erts a time-depende nt shear stress, arj_ = iuivY^'^ pv (Lan- 



dau fc Lifshitz 1959| ), that shakes the crust and can poten- 



tering, and / — 1/15. Finally, when both protons and 
neutrons are superfluid, electron-electron scattering yields 
/ = 5(p/pb). 

We begin by assuming that the crust is infinitely rigid, 
and hence the waves do not penetrate into the crust. In 
an oscillating flow, the thickness of the VBL, 5, is found 
by setting the oscillation frequency of the r-mode in the 
rotating frame, lo — 2mil/l{l + 1) = 217/3 {I — m = 2), to 
the inverse of the time it takes vorticity to diffuse across 
its width, i.e. S'^/u 



tially couple the crust to the core. If this coupling can drive 
the amplitude of the crustal motion to be comparable 
to the transverse motions in the r-modes ~ 3aR/2 for 
I — m — 2), then our picture of the core rubbing against a 
static crust would need to be reconsidered. The energy in 
an r-mode oscillation is ^ 10^^a^(i^s/l kHz)^ ergs, while 
the typical energy in torsional oscillations of the crust is 



roughly iC^a^ ergs ( [McDcrmott et al. 198^ ). Hence, if 



CO 



In the plane-parallel case. 
Landau fc Lifshitz (1959 ) show that the correction to the 
velocity in the VBL of an incompressible oscillatory flow 
falls off exponentially with a length 



LU 



1/2 



3 cm 



/1/2 fi kHz\ 

l^s J 



1/2 



(2) 



much smaller than R. This estimate neglects the Coriolis 
force. When the VBL is determined by the balance be- 
tween the Coriolis force, ^ 2riu, and viscosity, ~ vvjS^ 
(the Ekman problem), then the thickness is ^ ~ (t^/r2)^/^- 



the r-mode could couple efficiently to the crust, it would 
have adequate energy to drive crustal pulsations of compa- 
rable amplitude and potentially break the crust. We now 
estimate the magnitude of this coupling and show that, 
under most circumstances, it should be quite small. 

Let us look at a toy problem of a simple harmonic os- 
cillator (i.e., the crust) with a mode frequency loq and 
weak damping 7 <C wo driven by a harmonic force per 
unit mass /oe*"*. Away from resonances with crustal 
pulsation modes, the response of the crust is then just 
C± ~ /o/(wo — i^^)- Ignoring Coriolis force, the funda- 
mental I — 2 torsional mode has a frequency ~ 40 Hz, 
while the overtones have frequencies ^ 500n Hz, where 



roughly the same as eq. ( 2 ) . Hence, while the Coriolis force 
will change the angular dependence ot the velocity m the 
VBL, it is unlikely to change its thickness. 

The viscous dissipation can be handled quite simply for 
this thin VBL. Since all of the kinetic energy there is 
damped in one cycle, the Q of the oscillation is roughly 
just the ratio of the total volume to that in the VBL, 
or Q ~ RjU - Wi^T^ivsIl kUz)^/'^ / f^/'^ . The damping 
time is then a few hundred seconds, thus competitive with 
the growth time from gravitational wave emission, Tg„ — 
-146 s Mf 4i?6 ^(i^s/1 kHz)-^ where M1.4 = M/1.4M, 
and Re — R/10 km, of the / ~ % 
19981; iLindblom et al. 19981). 



n > 1 is the radial order of the overtone ( McDcrmott et al 



1988). Damping times o f these oscillations are on the or 



der of days to decades (McDcrmott et al. 1988), so the 
resonances are narrow. The spacing between the crustal 
modes, ~ 500 Hz, is comparable to the r-mode frequency. 



so on average, uj^ 



2 r-mode (Owen et al 



Q — w ~ w . Mass per unit area in the 
crust is Afcrust/47r_R^ = ph, where h ^ 1 km is the scale 
height at the base of the crust. Therefore, /o = <Jr±/ph 
and, in the likely situation that the r-mode frequency is not 
in resonance with a normal mode frequency of the crust, 
we find ^±/R ~ a{S/h) « IQ-^a/Tg. 

This simple model (and a preliminary calculation) shows 
that the induced toroidal displacements in the crust are 
much smaller than the r-mode displacements. In other 
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Fig. 1. — Critical spin frequencies for the r-mode instability. 
Solid line shows the critical frequency set by the viscous boundary 
layer and internal dissipation, for a star with normal nucleons in the 
core. The shading around the solid line displays the effect of core 
superfluidity. The dashed line shows the critical frequency as calcu- 
lated previously, neglecting the boundary layer friction. The LMXBs 
reside in the rectangular shaded region. The arrows show the fre- 
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Fig. 2. — Evolutionary scenarios for rapidly rotating neutron 
stars with crusts. The thick solid line shows the critical spin fre- 
quency for the r-mode instability, including boundary layer friction. 
The two thin solid loops show the evolution of a neutron star in an 
LMXB, accreting at 1O~*M0 yr"-*^, for two different values of the 
r-mode saturation amplitude a^, as marked on the plot. The dotted 
line shows the spin-down and cooling of a newborn rapidly rotating 
neutron star. 



words, even though the r-mode has plenty of energy to 
drive crustal pulsations, it cannot effectively shake the 
crust. Hence, neglecting the coupling should not invali- 
date our estimate of the damping time Trub- 

3. THE R-MODE INSTABILITY 

We now find the region in [vs, T) space where a NS is un- 
stable to the r-mode instability by balancing the viscous 
damping with the excitation due to gravitational radia- 
tion. When the crust is present, the r-mode is unstable 



for T,., b < |v 



The shaded box in Figure |] indicates the region of the 
{I's, T) space where the observed LMXBs (250 Hz < ^ 
600 Hz, T ~ (1 - 8) X IQS K; [Brown fc Bildsten 1998| ) 
reside. The arrows on the left-hand side of the plot show 
frequencies and the upper limits on the core temperatures 
for the two fastest MSPs. Clearly, the strength of the VBL 
dissipation means that r-modes are not excited in the ac- 
creting and colder [T ^ 10^ K) NSs. This result removes 
the discrepancy between the existence of 1.6 ms (625 Hz) 
millisecond p ulsars and the r-mode instability (Andersson 



> 1070 Hz 



4/11 
1.4 



_^1/11 



14/11 y2/ll 



2/11 



ct al. 1999b) and alleviates the disagreement between the 
observed and predicted quiescent luminosities of NS tran- 



(5) 



sients (Brown & Ushomirsky 199!;). 

We now re-examin e the therm ogravitational runaway 
scenario for LMXBs ( [Levin 1999D . The two closed loops 
in Figure show the evolution of the spin and temperature 



Howe ver, other viscous mechanisms^in the core also con- ^ j ^g accreting at IO-^Mq yr"! (i.e, same as in 



tributp to the damping. In Figure [l[ we show the critical 

'gw -l/'rub — 1/' 



spin fJuqueiicy whcic l/'/pw — l/'^■||b~l/'^vis^ — 0, wheiewg- 
use previou sly publishe d calculations of shear ( Lindblo"iij 
et al. 1998) and bulk (Lindblom ct al. 1999) viscous r- 



mode damping for the interior viscous damping time Tvisc ■ 
The solid line corresponds to the case where all nucleons 
are normal (/ ~ 1), while the shading around it repre- 
sents the range in frequencies when either neutrons or all 
nucleons are superfluid. For comparison, the dashed line 
shows the instability curve that neglects the effect of the 
viscous boundary layer. Our curves presume that the crust 
is present for all temperatures less than ^ 2 x 10^° K. For 
temperatures higher than that, the bulk viscosity is the 
dominant damping mechanism. 



1999), but including the VBL damping. The evolution is 
qualitatively the same. First, the star spins up until it 
reaches the instability line. The r-mode amplitude then 
grows until saturation (at amplitude Ofg), heating the star, 
and then spins it down due to gravitational wave emis- 
sion. The r-mode eventually stabilizes and the star cools 
down. For the chosen a^, the thermal runaway spins the 
star down to 200 — 400 Hz. The final frequency is not 
sensitive to the initial temperature of the star, but only to 
the assumed . The duration of the active r-mode phase 
is < yr, while the duration of the cool-down at constant 
spin is ^ 10^ yrs. The spinup from ~ 400 to ~ 1100 Hz 
takes roughly 5 x 10® yrs. 
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Therefore, in agreement with Levin (1999| ), we do not 



expect that any of the currently observed LMXBs are in 
the active r-mode phase. Moreover, this evolution can- 
i iot explain the clu stering of LMXB spins around 300 Hz 
dvan der Klis 1999D , since the duration of the cool-down 
at constant spin phase (horizontal leg in Figure ^) is much 
shorter than that of the spin- up phase (the vertical leg). 
Other mechani sms, such as mass quadrupole radiation 
from the crust ( Bildsten 1998| ) can likely explain such clus- 
tering. 

What about newborn NSs? NSs born in supernovae are 
thought to rotate near breakup, and have initial temper- 
atures Ti « 10^1 K. They quickly cool, roughly according 
to Tg = {t/Tc)-^/^ for T < Tj, where 1 yr is th e 

Urea cooling time at 10^ K ([Shapiro fc Teukolsky 19831 ). 
During the initial cooling stage, the r-mode instability line 
is as calculated previously (dashed line in Figure |l|), and 
the spin-dow n evolution proceeds according to the scenario 
described by Owen et al. (1998 ). However, a solid crust 
forms when the NS has cooled to 10^° K (the exact 

temperature depends on the composition), after which we 
expect the evolution to be significantly altered. 

The dotted line in Figure shows th e evolution of an d 
r of a newborn NS, computed as in Owen et al. (1998), 
with initial spin i/sa — 1320 Hz, but including the VBL 
and the heating from it. The star cools from T.^ = 10^^ K 
and enters the instability region. The r-mode then grows, 
saturates at Ofg = 1, and spins the star down. We find 
that the final spin frequency of the star is i^f = 415 Hz in - 
cluding the VBL, rather than 120 Hz ( pwen et al. 19981 ). 
The spindown is also of shorter duration, lasting only 10 s, 
rathe r than the ~ 1 yr of the original calculations of pwer] 
et al. 1(199^ ). For shorter Urea cooling time Tc, the final 
frequency is somewhat higher. The effect of the heating 
due to the VBL on the cooling time is negligible. In our 
simulation, the r-mode amplitude grows only by a factor 
of 3 from its initial value by the time the star cools to 
Tni, so we do not expect the presence of the r-mode to af- 
fect crust formation. In general, the rapid spindown phase 
begins when the Urea cooling time at the current temper- 
ature, T/{dT/dt) (X T~^, exceeds the gravitational wave 
spindown time, fl/(dfl/dt) oc aj^. Therefore, unless Tm 
is much less than 10^*^ K, the crust will form while the 
r-mode amplitude is still rather small. 

The strain amplitude is a very shallow function of Vs 



during the spin-down stage, he oc ly^"^ ( Owen et al. 199"8| ) 



However, most signal-to-noise {S/N) in a gravitational 
wave detector is accumulated at low v^, where the source 
spends most of its time. Hence, the total S/N depends 
sensitively on the low-frequen cy cutoff of the sp in-down 
evolution. Using the results of Owen et al. (1998 ), we find 
that raising the cutoff frequency to 415 Hz from 120 Hz 
reduces the S/N from 8 to ~ 2.5 for "enhanced" LIGO 
and a NS at a distance of 20 Mpc, where the event rate is 
expected to be a few per year. 



4. CONCLUSIONS AND FUTURE WORK 

Our initial work on the viscous boundary layer between 
the oscillating fluid in the core and the co-rotating crust 
shows that the dissipation there is very large, making it 
the predominant r-mode damping mechanism when the 
crust is present. The smallest spin frequency which allows 
the r-mode to operate is about 500 Hz, nearly a factor of 
five higher than previous estimates. This substantially re- 
duces the parameter space for the instability to operate, 
especially for older, colder NSs, such as those accreting in 
binaries and millisecond pulsars. It resolves the discrep- 
ancy between the theoretical understanding of the r-mode 
instability and the observations of millisecond pulsars and 
LMXBs. 

In our estimates we presumed that the boundary layer is 
laminar. Of course, turbulence in the VBL would increase 
the viscosity there and move the instability to higher rota- 
tion frequencies. We also showed that the r-modes induce 
some transverse motions in the crust, though these mo- 
tions appear to be negligibly small compared to the r-mode 
amplitude. The Coriolis force is likely to change the details 
of the eigenfunctions in the boundary layer, but it is un- 
likely to change its gross structure. Finally, based on the 
large conductivity contrast between the core and the crust, 
we have presumed that the VBL heating evenly spreads 
throughout the isothermal core. This may be an oversim- 
plification, and a more detailed time-dependent calcula- 
tion of the heat transport will yield the radial tempera- 
ture profile. This is especially important for the accreting 
systems, where the crust melting temperature s are lower 
because of th e smaller mean nuclear charges (Hacnsel & 
Zdunik 1990| ). 

Can a weak magnetic field substantially change our 
arguments? The Ohmic diffusion time across the VBL 
is many years, much longer than the oscillation pe- 
riod. Hence, any magnetic field protruding into the 
fiuid from the crust will be pulled and sheared in a 
nearly ideal MHD response. The field produced by 
the shear is SB/B ^ aR/6, and hence the magnetic 
field does not compete with viscous shear when B <C 
ip^n)^/^ « 10" G {us/l kHz)i/Vr8- The surface dipole 
field strengths in LMXBs and MSPs are lower than this 
value. If the same can be said about the B field at the 
base of the crust, then the viscous boundary layer we have 
discussed should apply. Of course, for younger and poten- 
tially much more magnetic neutron stars, the story will 
be modified, as the restoring force from the stretched field 
lines needs to be self-consistently included. 

We thank Andrew Gumming, Yuri Levin, and Lee Lind- 
blom for helpful conversations. L.B. is a Gottrell Scholar of 
the Research Gorporation and G.U. is a Fannie and John 
Hertz Foundation Fellow. This work was supported by 
the National Science Foundation through Grants AST97- 
31632 and PHY 94-07194. 
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